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ABSTRACT 

We report on spectro-imaging observations of the Herbig-Haro 2 outflow 
with the ISOCAM camera onboard the Infrared Space Observatory (ISO). The 
[Ne II] 12.81/im and [Ne III] 15.55/im lines are detected only towards the jet 
working surface (HH 2H), consistent with the high excitation of this knot in the 
optical range, while H 2 pure rotational emission is found all over the shocked 
region HH 2. The low energy transition S(2) traces warm gas (T ~ 400 K) 
peaked towards knots E-F and extended ejecta (T ~ 250 — 380 K) with masses of 
a few 10~ 3 M in the high-velocity CO outflow extending between the powering 
source and HH 2. Such emission could arise from low- velocity C-type shocks 
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(v ~ 10 — 15kms _1 ). The higher transitions S(3)-S(7) trace the emission of 
hot shocked gas (T = 1000 — 1400 K) from individual optical knots in the HH 2 
region. The ortho to para (OTP) ratio exhibits large spatial variations between 
1.2 (E) and 2.5 (H), well below its value at LTE. The emission of the S(3)- 
S(7) lines is well accounted for by planar C-shock models with a typical velocity 
V s = 20 — 30 km s _1 propagating into a medium of density ni = 10 4 — 10 5 cm" 3 
with an initial OTP ratio close to 1 in the pre-shock gas. In the leading edge of 
the jet, where the geometry of the emission allows a simple modelling, a good 
agreement is found with velocities derived from the optical proper motions mea- 
sured in the ionized gas. 

Subject headings: ISM: Herbig-Haro objects — ISM: individual (HH 1/2) — ISM: 
jets and outflows — ISM: molecules - stars: formation 



1. Introduction 

Bipolar outlows from embedded young stellar objects (YSOs) are perhaps one of the 
most spectacular manifestations of the star formation process. One of the best studied 
and brightest outflows is that of the Herbig-Haro (hereafter HH) 1/2 system (Reipurth 
1993). The HH 1/2 system lies in the Orion molecular cloud at 440 pc and subtends a ~ 2' 
angle. The VLA 1 embedded source (Pravdo et al. 1985) drives a highly collimated jet that 
reaches atomic gas velocities of ~ 480 km s" 1 (Eisloffel et al. 1994) and produces shocks 
of ~ 100 km s" 1 at its main working surface or bow shock ( Noriega- Crespo et al. 1989). At 
optical wavelengths there are at least 3 distinct jet flows arising within 5" of the VLA 1 source 
and the HH 1/2 bow shocks display a complex morphology in high spatial resolution HST 
images (Bally et al. 2002). The optical jet is associated with a molecular outflow whose high- 
velocity (the "molecular jet") covers deprojected velocities of 15 — 80kms _1 with respect to 
the ambient cloud (Moro-Martm et al. 1999). Because of these characteristics the HH 1/2 
system is a perfect target to study the gas properties in a shock heated environment as 
well as the spatial distribution and the nature of these shocks. One of the promises of ISO 
was precisely to be able to discern between the different shocks, either C-type or J-type, 
occurring in a molecular environment. 



1 Bascd on observations with ISO, an ESA project with instruments funded by ESA Member States 
(especially the PI countries: France Germany, the Netherlands and the United Kingdom) and with the 
participation of ISAS and NASA 
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We present here spectro-imaging observations of the VLA 1 counter-jet and the HH 2 
region obtained between 5 and 17/im with the ISOCAM camera onboard ISO. We find that 
the H 2 pure rotational lines S(2) to S(7) arise from two physically distinct regions : a faint 
extended warm gas component associated with the high-velocity CO jet, seen mainly in S(2), 
and several more compact peaks tracing hotter shocked gas in the individual knots A-L of 
HH 2 detected in all lines. We report also on the presence of the [Ne III] 15.55/zm and 
[Ne II] 12.8/im lines. HH 2 is the only outflow of low-luminosity where these infrared lines 
have been detected so far (see Cabrit et al. 1998). 



2. Observations 

All observations were obtained with the ISO satellite (Kessler et al., 1996) and the 
ISOCAM instrument (Cesarsky et al. 1996). The low resolution spectra (A/AA = 40) 
between 5 and 17/im were obtained in revolution 691 with the Circular Variable Filter (CVF) 
with a pixel scale of 6" and a total field of view of 3' centered on the HH 2 object. The last 
pipeline version of the data (OLP10) has been processed following the package developed 
at the Institut d'Astrophysique Spatiale, which removes reasonably well the problem of 
transients. The size (HPFW) of the Point Spread Function (PSF) is 6" for a pixel scale 
of 6". 

In order to establish accurate astrometry, we used a second CVF map containing the 
optically visible Cohen-Schwartz (CS) star, taken in revolution 873 with a 3" pixel scale. 
This second dataset was processed in the same way as mentioned above. Unfortunately, 
internal reflections between the CVF and the field lens produced spurious ghosts of the CS 
star, which is a strong IR emitter, over the full field, preventing any quantitative analysis of 
the H 2 line emission in the second data set. However, the presence of the VLA 1 protostar 
in both data cubes allows to derive an accurate astrometry for the first CVF image taken 
with a 6" pixel scale. The data are presented in Figure 1. Coordinates are offsets (arcsec) 
relative to the position of VLAl : a 200 o = 05 h 36 m 22.6 s , 5 2 ooo = -06°46'25". The interstellar 
extinction towards HH 2 was estimated by Hartmann & Raymond (1984), who measured 
typical reddenings E(B — V) = 0.11 — 0.44. Based on the extinction curve of Rieke & 
Lebofsky (1985), it appears that the flux dereddening corrections are negligible and we use 
uncorrected flux values in what follows. 
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3. Results 

Individual CVF spectra towards various optical knots in HH 2 are displayed in the 
panels of Fig. lb. Knot positions are referred by letters A to L following the nomenclature 
of Eisloeffel et al. (1994). At most positions, the bulk of the emission in the 5-17/zm range 
comes from the H 2 pure rotational lines S(2) to S(7). An exception is the region near HH 2 H 
(positions H,B,D) where the fine structure ionic lines [Ne II] 12.81/wn and [Ne III] 15.55 fxm 
are also detected, and even dominate over H 2 lines at the nominal position of knot H. The 
spatial distribution of emission flux in the H 2 v=0-0 S(2), S(3), S(5) lines and in the [Ne 
II] 12.8 fim ionic line is illustrated in Fig. la and Fig. lc. Three types of morphologies are 
observed, depending on the line excitation : 

(i) The ionic lines [Ne II] and [Ne III] show a single peak towards knot H (of typical 
size ~ 7" at HPFW). (ii) The intermediate excitation lines S(3) and S(5) show two peaks 
of comparable brightness, one encompassing E, and another peak centered between knots H 
and D, shifted by 3" from the ionic line peak, (iii) The low energy S(2) line in Fig. la shows 
a single peak (of size 13") encompassing knots E and F. The emission peaks at the tip of the 
molecular jet 15" downstream the CO brightness peak. The lowest contours of S(2) emission 
reveal a broad pedestal that points towards VLA1/4 (Fig. la). Note that the first contour is 
at 5a above the map noise level. The pedestal overlaps well with the CO "jet" and suggests 
that both the CO and H 2 emissions are somewhat related, although tracing different regions 
in the "jet". 

The overall agreement between the brightness distributions of the pure H 2 rotational 
lines and the higher excitation line 2.12/zm 1-0 S(l) (Davis et al. 1994) is good. Comparison 
of the line intensities indicates that the H 2 pure rotational lines are collisionally excited 
(Wolfire & Konigl 1991). We show in the right panels of Figure 1 the H 2 rotational diagrams 
obtained towards three representative positions : a region in the CO jet (polygon in Fig. la), 
knot HH 2E (peak of H 2 no ionic line emission) and knot HH 2D (both ionic and molecular 
features). Overall, the fluxes of the lines S(J > 3) match well a linear fit at all positions. 
The excitation temperature T cx , the ortho to para ratio (OTP), and the total H 2 column 
density of this "hot" component are estimated from a common linear fit to the temperature 
for the ortho and para species in the excitation diagram, and are summarized in Table 1. 

Our analysis shows evidence for a hot gas layer with a temperature in the range 970- 
1400 K, total H 2 column density ranging between 0.18 (H) and 1.6 x 10 19 cm -2 (E), and 
OTP ratio ranging between 1.2 (E) and 2.5 (H). There is no spatial trend in the values, i.e. 
upstream knots do not have a higher OTP ratio. Interestingly, there appears to be some 
correlation between variations in OTP and T cx : the smallest OTP (1.2) is found towards 
HH 2E, which has among the lowest T ex , while the highest OTP (2.5) is found towards HH 
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2H, which has the highest T cx . Knots A,C,D, and L fall in between these two extremes. 
Knots B,K,G are the only positions that deviate from this trend (low T cx but OTP close to 
2). The OTP ratios measured are systematically lower than the LTE value at all positions 
but H. It is consistent with the low-excitation temperatures measured in the hot gas if H 2 is 
excited by means of a C-shock (see e.g. Wilgenbus et al. 2000). 

The rotational diagrams show that, at each position, the population of the upper level of 
the S(2) line lies well above the hot component fit to the other transitions. This implies that 
the pedestal H 2 emission is dominated by a second gas component at a lower temperature, 
discussed in Sect. 4.3 



4. Discussion 

4.1. Ionic emission in HH 2H 

The detection of the [Ne II] 12.8 /im line towards HH 2H is a clear signature of J-shocks 
with velocities above GOkms^ 1 (Hollenbach & McKee 1989). The presence of [Ne III] is 
consistent also with the detection of the [Si II] 34.8 fxm line in HH 2, which suggests shock 
velocities of 100 - 140 km s* 1 (Molinari & Noriega-Crespo 2002). That [Ne II] and [Ne III] 
are detected only towards HH 2H is in line with the particularly high excitation of this knot 
in the optical (Bohm & Solf 1992), and its proper motion larger than 400 km s^ 1 (Bally et 
al. 2002). The high excitation of HH 2H is attributed to it being the current location of the 
jet working surface (Bally et al. 2002). 

The shock speeds inferred in knot H exceed by far the speed at which molecules are 
dissociated (~ 70kms _1 for a C-shock at n (H) = 10 4 cm~ 3 ; LeBourlot et al. 2002, and 
~ 25kms~ 1 for a J-shock; Hollenbach & McKee 1989). This is consistent with the H 2 peak 
in S(5) being spatially shifted from HH 2H. A similar shift was seen in the H 2 1-0 S(l) line 
(Noriega-Crespo & Garnavich, 1994). It indicates that the H 2 comes from a separate, lower 
velocity shock not physically associated with knot H. 

4.2. Shock-excitation of the hot H 2 component in HH knots 

Models of non-dissociative J-shocks (Wilgenbus et al. 2000) or dissociative J-shocks 
with H 2 reformation (Flower et al. 2003) fail to reproduce the CVF data as the predicted 
line intensities, expecially S(3)-S(5), are much too weak, or the required shock velocity is 
too low (~ lOkms -1 ). On the contrary, everywhere in HH 2, the recent models of planar 
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C-shocks (Le Bourlot et al. 2002; Cabrit et al. 2003) provide a much better match with the 
observations. 

Comparison with the observations favors models with velocity in the range 20—30 km s" 1 
and a preshock density in the range 10 5 — 10 4 cm -3 respectively (see the rotational diagrams 
in Fig. 1). Observations of higher- J H 2 lines would allow to better constrain the density in 
the pre-shock gas. This range of densities also agrees with the determinations obtained in 
the ambient molecular cloud (Girart et al. 2002). The OTP ratio in the pre-shock gas is 
found ~ 1. The observed variations in the OTP ratio from knot to knot can be produced 
by changes in the shock speed or preshock density, and do not require variations in the 
initial OTP ratio in the preshock gas. The characteristic size of the shock emitting region 
is constrained to be of the order of 0.8" — 1.7" depending on the individual knots, which 
compares well with their size in the optical. 

The above shock velocities are derived from planal models whereas HST images provide 
large evidence for numerous bow-shaped features in HH 2 (Bally et al. 2002). For several 
knots near the leading edge of HH 2 where the geometry is rather simple and bow-shaped 
features easily identified (e.g. F, E, L), the shock velocities compare well with the proper 
motions of the bow wings in the optical (e.g. 50 — 60 km s" 1 at E) taking into account the 
obliquity of the shock with respect to the direction of the motion. 

The absence of spatial trend in the OTP values probably results from the various shocks 
associated with the mini bows, unresolved in our CVF images. Each of these shocks modify 
the OTP ratio depending on its own excitation conditions. 

4.3. Warm extended H 2 component 

An upper limit to the temperature of the extended warm component in the H 2 pedestal 
can be derived from the ratio of the S(2) and S(3) lines (adopting an OTP ratio similar to 
that in the ambient cloud, close to 1). We find T < 420 K at the H 2 peak and T < 380 K 
towards the CO jet. Note that a similar constraint, independent of the OTP ratio, is obtained 
by assuming that at most 50% of the S(4) flux comes from the pedestal. 

Conversely, in the optically thin limit, we can also derive a lower limit to the temperature 
by imposing that the warm component contributes to most of the observed flux of the J= 
2-1 line in the CO "jet" and assuming a standard abundance [CO]/[H 2 ] = 10~ 4 (Flower and 
Pineau des Forets 1994 showed indeed that the CO abundance varies little in C-shocks). 
This temperature is found by equating expressions of the H 2 column density derived from 
both tracers. Towards the region of maximum CO flux (10 — 20 K km s" 1 ) in the jet the H 2 
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S(2) line flux is ~ 8 mJy/pixel, which yields a minimum temperature of ~ 240 — 300 K in 
the pedestal. This determination depends very little on the adopted OTP value. At the H 2 
peak the CO(2-l) emission is much weaker (~ 6Kkms~ 1 ) and we infer T > 400 K. 

The warm gas at 240 — 420 K detected in the pure H 2 rotational lines of low-J is 
predominantly concentrated towards the tip of the high-velocity CO jet. The S(2) line 
suggests a warm gas column density of about 2 x 10 20 cm -2 at knots E-F and a corresponding 
mass of ~ 2 x 10~ 3 M . The hot gas seen towards the high- velocity CO jet has a 10 times 
lower column density (Fig. 1). To reproduce the observed J=4 population, as well as the 
limits on the excitation temperature, a low-velocity C-shock with a filling factor of 1, implied 
by the extended character of the S(2) emission, is needed. Then, a rather slow shock at 
a velocity of ~ 10 — lSkms" 1 into gas of density n(H) ~ 10 5 cm -3 , with an excitation 
temperature of 250 K would account for the H 2 emission of the pedestal along the CO jet. 

Observations of higher J CO rotational lines at higher angular resolution are needed to 
determine more precisely the physical conditions of the high-velocity CO gas and to compare 
it with that traced in S(2) and with molecular shock predictions. 

J. Cernicharo acknowledges Spanish DGES for this research under grants PNAYA 2000- 
1784 and ESP2001-4516. ANC's research was carried out at the Jet Propulsion Laboratory, 
California Institute of Technology, under a contract with NASA; and partially supported by 
NASA-APD Grant NRA0001-ADP-096. 
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Table I. Parameters of the hot molecular gas derived from rotational diagram analysis. 
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Fig. 1 - (left) (a) : Intensity contour map of the H2 0-0 S(2) line (pink) superposed on a [S II] 
color image of the HH 1/2 region (Reipurth et al. 1993). Contour levels are 4, 7 10, 15, 20, 
45mJy/pixel. We have superposed the emission of the high-velocity CO outflow in white con- 
tours. First contour and contour interval are 5 and 2.5 K km s _1 respectively. The position of the 
VLA 1-4 sources is indicated by yellow filled circles. Black dots mark the location of the knots A-L. 
(bl-b9) : CVF spectra at selected positions. Fluxes are in mJy/pixel. The jet spectrum is an 
average of the individual CV spectra between HH 2 and VLA 1 in the area of the polygon (yellow) 
on panel a. The wavelength of the pure rotational H 2 transitions S(2) to S(7) is marked with red 
ticks, (c) : Superposed on the same [S II] picture : emission contour map of the lines H2 S(3) (cl), 
H 2 S(5) (c2) and [Ne II] (c3). Contour levels are 0.1, 0.2, ... 0.9 times the brightness peak. The H 2 
intensity contour maps in panels (a) and (c) have been convolved with a gaussian of 9" HPFW. 
(right) Rotational diagrams for the H 2 rotational states J=3-9. The logarithm of Nj/(gjg s ) is plot- 
ted against Ej/kb, where Nj is the beam-averaged column density, gj is the rotational degeneracy, 
g s is the spin degeneracy, and Ej is the energy of the rotational level J. Open squares mark the 
data points. Blue straight lines show the best fit for a gas component of uniform temperature. 
The gas properties are given in each panel. The dashed lines in each panel correspond to the same 
temperature although they differ in the total column density for the orto and para species. Red 
dashed lines trace the fit to the warm extended component (pedestal). Green filled stars show the 
best planar C-shock model that accounts for the J > 3 data points. The parameters of the model 
(shock velocity V s , hydrogen nuclei density no(H) and ortho to para ratio otpo in the pre-shock 
gas) and the pixel filling factor ff are given for HH 2E and HH 2D. Below the panels of HH 2E and 
HH 2D are given the best fit parameters for planar C-shock models. 
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ABSTRACT 



We report on spectro-imaging observations of the Herbig-Haro 2 outflow 
with the ISOCAM camera onboard the Infrared Space Observatory (ISO). The 
[Ne II] 12.81/ira and [Ne III] 15.55/ira lines are detected only towards the jet 
working surface (HH 2H), consistent with the high excitation of this knot in the 
optical range, while H 2 pure rotational emission is found all over the shocked 
region HH 2. The low energy transition S(2) traces warm gas (T ~ 400 K) 
peaked towards knots E-F and extended ejecta (T ~ 250 — 380 K) with masses of 
a few 10~ 3 M & in the high-velocity CO outflow extending between the powering 
source and HH 2. Such emission could arise from low- velocity C-type shocks 
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(v ~ 10 — 15kms _1 ). The higher transitions S(3)-S(7) trace the emission of 
hot shocked gas (T = 1000 — 1400 K) from individual optical knots in the HH 2 
region. The ortho to para (OTP) ratio exhibits large spatial variations between 
1.2 (E) and 2.5 (H), well below its value at LTE. The emission of the S(3)- 
S(7) lines is well accounted for by planar C-shock models with a typical velocity 
V s = 20 — 30kms _1 propagating into a medium of density n; = 10 4 — 10 5 cm -3 
with an initial OTP ratio close to 1 in the pre-shock gas. In the leading edge of 
the jet, where the geometry of the emission allows a simple modelling, a good 
agreement is found with velocities derived from the optical proper motions mea- 
sured in the ionized gas. 

Subject headings: ISM: Herbig-Haro objects — ISM: individual (HH 1/2) — ISM: 
jets and outflows — ISM: molecules - stars: formation 



1. Introduction 

Bipolar outlows from embedded young stellar objects (YSOs) are perhaps one of the 
most spectacular manifestations of the star formation process. One of the best studied 
and brightest outflows is that of the Herbig-Haro (hereafter HH) 1/2 system (Reipurth 
1993). The HH 1/2 system lies in the Orion molecular cloud at 440 pc and subtends a ~ 2' 
angle. The VLA 1 embedded source (Pravdo et al. 1985) drives a highly collimated jet that 
reaches atomic gas velocities of ~ 480 km s -1 (Eisloffel et al. 1994) and produces shocks 
of ~ 100 km s -1 at its main working surface or bow shock (Noriega-Crespo et al. 1989). At 
optical wavelengths there are at least 3 distinct jet flows arising within 5" of the VLA 1 source 
and the HH 1/2 bow shocks display a complex morphology in high spatial resolution HST 
images (Bally et al. 2002). The optical jet is associated with a molecular outflow whose high- 
velocity (the "molecular jet") covers deprojected velocities of 15 — 80 kms -1 with respect to 
the ambient cloud (Moro-Martin et al. 1999). Because of these characteristics the HH 1/2 
system is a perfect target to study the gas properties in a shock heated environment as 
well as the spatial distribution and the nature of these shocks. One of the promises of ISO 
was precisely to be able to discern between the different shocks, either C-type or J-type, 
occurring in a molecular environment. 



1 Based on observations with ISO, an ESA project with instruments funded by ESA Member States 
(especially the PI countries: France Germany, the Netherlands and the United Kingdom) and with the 
participation of ISAS and NASA 
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We present here spectro-imaging observations of the VLA 1 counter-jet and the HH 2 
region obtained between 5 and 17/im with the ISOCAM camera onboard ISO. We find that 
the H 2 pure rotational lines S(2) to S(7) arise from two physically distinct regions : a faint 
extended warm gas component associated with the high-velocity CO jet, seen mainly in S(2), 
and several more compact peaks tracing hotter shocked gas in the individual knots A-L of 
HH 2 detected in all lines. We report also on the presence of the [Ne III] 15.55/im and 
[Ne II] 12.8/ira lines. HH 2 is the only outflow of low-luminosity where these infrared lines 
have been detected so far (see Cabrit et al. 1998). 



2. Observations 

All observations were obtained with the ISO satellite (Kessler et al., 1996) and the 
ISOCAM instrument (Cesarsky et al. 1996). The low resolution spectra (A/AA = 40) 
between 5 and 17/im were obtained in revolution 691 with the Circular Variable Filter (CVF) 
with a pixel scale of 6" and a total field of view of 3' centered on the HH 2 object. The last 
pipeline version of the data (OLP10) has been processed following the package developed 
at the Institut d'Astrophysique Spatiale, which removes reasonably well the problem of 
transients. The size (HPFW) of the Point Spread Function (PSF) is 6" for a pixel scale 
of 6". 

In order to establish accurate astrometry, we used a second CVF map containing the 
optically visible Cohen-Schwartz (CS) star, taken in revolution 873 with a 3" pixel scale. 
This second dataset was processed in the same way as mentioned above. Unfortunately, 
internal reflections between the CVF and the field lens produced spurious ghosts of the CS 
star, which is a strong IR emitter, over the full field, preventing any quantitative analysis of 
the H 2 line emission in the second data set. However, the presence of the VLA 1 protostar 
in both data cubes allows to derive an accurate astrometry for the first CVF image taken 
with a 6" pixel scale. The data are presented in Figure 1. Coordinates are offsets (arcsec) 
relative to the position of VLA1 : a 200 o = 05 /l 36 m 22.6 s , 8 20 oo = -06°46'25". The interstellar 
extinction towards HH 2 was estimated by Hartmann & Raymond (1984), who measured 
typical reddenings E(B — V) = 0.11 — 0.44. Based on the extinction curve of Rieke & 
Lebofsky (1985), it appears that the flux dereddening corrections are negligible and we use 
uncorrected flux values in what follows. 
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3. Results 

Individual CVF spectra towards various optical knots in HH 2 are displayed in the 
panels of Fig. lb. Knot positions are referred by letters A to L following the nomenclature 
of Eisloeffel et al. (1994). At most positions, the bulk of the emission in the 5-17 fira range 
comes from the H 2 pure rotational lines S(2) to S(7). An exception is the region near HH 2 H 
(positions H,B,D) where the fine structure ionic lines [Ne II] 12.81/ira and [Ne III] 15.55 fira 
are also detected, and even dominate over H 2 lines at the nominal position of knot H. The 
spatial distribution of emission flux in the H 2 v=0-0 S(2), S(3), S(5) lines and in the [Ne 
II] 12.8 fira ionic line is illustrated in Fig. la and Fig. lc. Three types of morphologies are 
observed, depending on the line excitation : 

(i) The ionic lines [Ne II] and [Ne III] show a single peak towards knot H (of typical 
size 7" at HPFW). (ii) The intermediate excitation lines S(3) and S(5) show two peaks 
of comparable brightness, one encompassing E, and another peak centered between knots H 
and D, shifted by 3" from the ionic line peak, (iii) The low energy S(2) line in Fig. la shows 
a single peak (of size 13") encompassing knots E and F. The emission peaks at the tip of the 
molecular jet 15" downstream the CO brightness peak. The lowest contours of S(2) emission 
reveal a broad pedestal that points towards VLA1/4 (Fig. la). Note that the first contour is 
at 5<T above the map noise level. The pedestal overlaps well with the CO "jet" and suggests 
that both the CO and H 2 emissions are somewhat related, although tracing different regions 
in the "jet". 

The overall agreement between the brightness distributions of the pure H 2 rotational 
lines and the higher excitation line 2.12/ira 1-0 S(l) (Davis et al. 1994) is good. Comparison 
of the line intensities indicates that the H 2 pure rotational lines are collisionally excited 
(Wolfire & Konigl 1991). We show in the right panels of Figure 1 the H 2 rotational diagrams 
obtained towards three representative positions : a region in the CO jet (polygon in Fig. la), 
knot HH 2E (peak of H 2 no ionic line emission) and knot HH 2D (both ionic and molecular 
features). Overall, the fluxes of the lines S(J > 3) match well a linear fit at all positions. 
The excitation temperature T ex , the ortho to para ratio (OTP), and the total H 2 column 
density of this "hot" component are estimated from a common linear fit to the temperature 
for the ortho and para species in the excitation diagram, and are summarized in Table 1. 

Our analysis shows evidence for a hot gas layer with a temperature in the range 970- 
1400 K, total H 2 column density ranging between 0.18 (H) and 1.6 X 10 19 cm -2 (E), and 
OTP ratio ranging between 1.2 (E) and 2.5 (H). There is no spatial trend in the values, i.e. 
upstream knots do not have a higher OTP ratio. Interestingly, there appears to be some 
correlation between variations in OTP and T ex : the smallest OTP (1.2) is found towards 
HH 2E, which has among the lowest T ex , while the highest OTP (2.5) is found towards HH 
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2H, which has the highest T ex . Knots A,C,D, and L fall in between these two extremes. 
Knots B,K,G are the only positions that deviate from this trend (low T ex but OTP close to 
2). The OTP ratios measured are systematically lower than the LTE value at all positions 
but H. It is consistent with the low-excitation temperatures measured in the hot gas if H 2 is 
excited by means of a C-shock (see e.g. Wilgenbus et al. 2000). 

The rotational diagrams show that, at each position, the population of the upper level of 
the S(2) line lies well above the hot component fit to the other transitions. This implies that 
the pedestal H 2 emission is dominated by a second gas component at a lower temperature, 
discussed in Sect. 4.3 



4. Discussion 

4.1. Ionic emission in HH 2H 

The detection of the [Ne II] 12.8 jiva line towards HH 2H is a clear signature of J-shocks 
with velocities above 60kms _1 (Hollenbach & McKee 1989). The presence of [Ne III] is 
consistent also with the detection of the [Si II] 34.8 jiva line in HH 2, which suggests shock 
velocities of 100 - 140 km s" 1 (Molinari k Noriega-Crespo 2002). That [Ne II] and [Ne III] 
are detected only towards HH 2H is in line with the particularly high excitation of this knot 
in the optical (Bohm & Solf 1992), and its proper motion larger than 400 km s -1 (Bally et 
al. 2002). The high excitation of HH 2H is attributed to it being the current location of the 
jet working surface (Bally et al. 2002). 

The shock speeds inferred in knot H exceed by far the speed at which molecules are 
dissociated (~ 70kms _1 for a C-shock at n (H) = 10 4 cm -3 ; LeBourlot et al. 2002, and 
~ 25kms _1 for a J-shock; Hollenbach & McKee 1989). This is consistent with the H 2 peak 
in S(5) being spatially shifted from HH 2H. A similar shift was seen in the H 2 1-0 S(l) line 
(Noriega-Crespo & Garnavich, 1994). It indicates that the H 2 comes from a separate, lower 
velocity shock not physically associated with knot H. 

4.2. Shock-excitation of the hot H 2 component in HH knots 

Models of non-dissociative J-shocks (Wilgenbus et al. 2000) or dissociative J-shocks 
with H 2 reformation (Flower et al. 2003) fail to reproduce the CVF data as the predicted 
line intensities, expecially S(3)-S(5), are much too weak, or the required shock velocity is 
too low (~ lOkms -1 ). On the contrary, everywhere in HH 2, the recent models of planar 
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C-shocks (Le Bourlot et al. 2002; Cabrit et al. 2003) provide a much better match with the 
observations. 

Comparison with the observations favors models with velocity in the range 20 — 30 km s -1 
and a preshock density in the range 10 5 — 10 4 cm -3 respectively (see the rotational diagrams 
in Fig. 1). Observations of higher- J H 2 lines would allow to better constrain the density in 
the pre-shock gas. This range of densities also agrees with the determinations obtained in 
the ambient molecular cloud (Girart et al. 2002). The OTP ratio in the pre-shock gas is 
found ~ 1. The observed variations in the OTP ratio from knot to knot can be produced 
by changes in the shock speed or preshock density, and do not require variations in the 
initial OTP ratio in the preshock gas. The characteristic size of the shock emitting region 
is constrained to be of the order of 0.8" — 1.7" depending on the individual knots, which 
compares well with their size in the optical. 

The above shock velocities are derived from planal models whereas HST images provide 
large evidence for numerous bow-shaped features in HH 2 (Bally et al. 2002). For several 
knots near the leading edge of HH 2 where the geometry is rather simple and bow-shaped 
features easily identified (e.g. F, E, L), the shock velocities compare well with the proper 
motions of the bow wings in the optical (e.g. 50 — 60kms _1 at E) taking into account the 
obliquity of the shock with respect to the direction of the motion. 

The absence of spatial trend in the OTP values probably results from the various shocks 
associated with the mini bows, unresolved in our CVF images. Each of these shocks modify 
the OTP ratio depending on its own excitation conditions. 

4.3. Warm extended H 2 component 

An upper limit to the temperature of the extended warm component in the H 2 pedestal 
can be derived from the ratio of the S(2) and S(3) lines (adopting an OTP ratio similar to 
that in the ambient cloud, close to 1). We find T < 420 K at the H 2 peak and T < 380 K 
towards the CO jet. Note that a similar constraint, independent of the OTP ratio, is obtained 
by assuming that at most 50% of the S(4) flux comes from the pedestal. 

Conversely, in the optically thin limit, we can also derive a lower limit to the temperature 
by imposing that the warm component contributes to most of the observed flux of the J= 
2-1 line in the CO "jet" and assuming a standard abundance [CO]/[H 2 ] = 10~ 4 (Flower and 
Pineau des Forets 1994 showed indeed that the CO abundance varies little in C-shocks). 
This temperature is found by equating expressions of the H 2 column density derived from 
both tracers. Towards the region of maximum CO flux (10 — 20 K kms -1 ) in the jet the H 2 
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S(2) line flux is ~ 8mJy/pixel, which yields a minimum temperature of ~ 240 — 300 K in 
the pedestal. This determination depends very little on the adopted OTP value. At the H 2 
peak the CO(2-l) emission is much weaker (~ 6 Kkms -1 ) and we infer T > 400 K. 

The warm gas at 240 — 420 K detected in the pure H 2 rotational lines of low-J is 
predominantly concentrated towards the tip of the high-velocity CO jet. The S(2) line 
suggests a warm gas column density of about 2 X 10 2 ° cm -2 at knots E-F and a corresponding 
mass of ~ 2 X 10~ 3 M & . The hot gas seen towards the high- velocity CO jet has a 10 times 
lower column density (Fig. 1). To reproduce the observed J=4 population, as well as the 
limits on the excitation temperature, a low- velocity C-shock with a filling factor of 1, implied 
by the extended character of the S(2) emission, is needed. Then, a rather slow shock at 
a velocity of ~ 10 — 15kms _1 into gas of density n(H) ~ 10 5 cm -3 , with an excitation 
temperature of 250 K would account for the H 2 emission of the pedestal along the CO jet. 

Observations of higher J CO rotational lines at higher angular resolution are needed to 
determine more precisely the physical conditions of the high- velocity CO gas and to compare 
it with that traced in S(2) and with molecular shock predictions. 
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1784 and ESP2001-4516. ANC's research was carried out at the Jet Propulsion Laboratory, 
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Table 1. Parameters of the hot molecular gas derived from rotational diagram analysis. 
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Fig. 1 - (left) (a) : Intensity contour map of the H2 0-0 S(2) line (pink) superposed on a [S II] 
color image of the HH 1/2 region (Reipurth et al. 1993). Contour levels are 4, 7 10, 15, 20, 
45mJy/pixel. We have superposed the emission of the high-velocity CO outflow in white con- 
tours. First contour and contour interval are 5 and 2.5 K km s _1 respectively. The position of the 
VLA 1-4 sources is indicated by yellow filled circles. Black dots mark the location of the knots A-L. 
(bl-b9) : CVF spectra at selected positions. Fluxes are in mJy/pixel. The jet spectrum is an 
average of the individual CV spectra between HH 2 and VLA 1 in the area of the polygon (yellow) 
on panel a. The wavelength of the pure rotational H2 transitions S(2) to S(7) is marked with red 
ticks, (c) : Superposed on the same [S II] picture : emission contour map of the lines H2 S(3) (cl), 
H2 S(5) (c2) and [Ne II] (c3). Contour levels are 0.1, 0.2, ... 0.9 times the brightness peak. The H2 
intensity contour maps in panels (a) and (c) have been convolved with a gaussian of 9" HPFW. 
(right) Rotational diagrams for the H2 rotational states J=3-9. The logarithm of Nj/(gjg s ) is plot- 
ted against Ej/kb, where Nj is the beam-averaged column density, gj is the rotational degeneracy, 
g s is the spin degeneracy, and Ej is the energy of the rotational level J. Open squares mark the 
data points. Blue straight lines show the best fit for a gas component of uniform temperature. 
The gas properties are given in each panel. The dashed lines in each panel correspond to the same 
temperature although they differ in the total column density for the orto and para species. Red 
dashed lines trace the fit to the warm extended component (pedestal). Green filled stars show the 
best planar C-shock model that accounts for the J > 3 data points. The parameters of the model 
(shock velocity V s , hydrogen nuclei density no(H) and ortho to para ratio otpo in the pre-shock 
gas) and the pixel filling factor ff are given for HH 2E and HH 2D. Below the panels of HH 2E and 
HH 2D are given the best fit parameters for planar C-shock models. 
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Fig. 1.- 
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